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Abstract 

We consider an extension of the MSSM with an added vectorlike top partner. Our aim is 
to revisit to what extent such an extension can raise the Higgs boson mass through radiative 
corrections and help ameliorate the MSSM hierarchy problem, and to specify what experimental 
probes at the LHC will find or exclude this possibility during the high-luminosity phase. Direct 
detection, precision electroweak and precision Higgs analyses are all commissioned to this end. To 
achieve the 125 GeV Higgs boson mass, we find that superpartner masses can be reduced by a 
factor of more than three in this scenario compared to the MSSM without the extra vectorlike top 
quark, and that during the high-luminosity phase of the LHC precision Higgs analysis is expected 
to become the most powerful experimental probe of the scenario. 
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I. INTRODUCTION 

Despite the lack of experimental confirmation from the hrst run of LHC, supersymmetry 
is still a promising solution to the Standard Model (SM) hierarchy problem. In the minimal 
realization of supersymmetry, the Higgs boson mass at tree level is bounded by the Z boson 
mass and needs to be lifted up by radiative corrections from superpartners. This calls for 
large superpartner masses that introduce a new hierarchy between the weak scale and the 
scale of supersymmetry. This is often called the little hierarchy problem of the MSSM. 

We will focus on an extension of MSSM with a vectorlike top quark partner. This is the 
simplest of vectorlike matter extensions [MQ] that can effectively reduce the little hierarchy 
due to large new contributions it induces to the Higgs mass. To illustrate this point we use 
the simplest possible supersymmetry spectrum with all soft terms at the scale Msusy- The 
only exceptions are the A terms equal to —Msusy- Also, the Higgs boson soft masses and 
B parameters are chosen to accommodate correct electroweak symmetry breaking. To this 
very simple spectrum we add a vectorlike top multiplet, t' and F, where t' has the quantum 
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numbers of the right-handed top quark and i' is its conjugate. The soft masses of the 
scalar components of t' and i' are also equal to Msusy- We include these new superhelds 
into the superpotential and calculate the contribution to the Higgs boson mass. As shown 
previously in different contexts as well m, the addition of vectorlike states that mix with 
the MSSM fermions can raise the Higgs boson mass, thereby enabling smaller superpartner 
masses to achieve nih = 125 GeV through these additional radiative corrections. 

We then calculate possible experimental exclusions or detections coming from precision 
electroweak measurements, corrections to Higgs boson properties and direct detection of the 
new vectorlike state. Finally we compare the impact of all these bounds on our model and 
calculate the lowest possible Msusy consistent with these bounds. 

One key result is that the most constraining of the three experimental analyses is usually 
the modihcation of Higgs boson properties, except when there is large tan/? and small 
mixing. In that case, the direct detection of exotic vectorlike states at the LHC can be 
more probing. We also will show that when including all the constraints, Msusy can still be 
lowered 3 to 5 times compared to the MSSM and still yield rrih = 125 GeV. Thus even a very 
simple vectorlike quark extension can signihcantly ameliorate the little hierarchy problem of 
MSSM. Its important to point out that since we do not consider a specihc UV completion, 
our measure of the little hierarchy problem is simply the splitting between the electroweak 
scale and the SUSY scale, rather than a result of some specihc hne-tuning measure. 


II. MSSM WITH VECTORLIKE TOP PARTNER 

The superpotential of the MSSM with an additional vectorlike top partner (omitting 
small Yukawa couplings of the hrst two families), reads 


W = YtQHj + Yt,QHj + mt't + Mt't' + YbQHj) + Y^LHdf + 


( 1 ) 


The above superpotential leads to the following mass matrix in the basis T = 


M. = 


0 mt 

m| 0 


nii = 


YtV2 Yt'V2 
m M 


( 2 ) 


where v = ~ 246, tan/? = V 2 /V 1 and V 2 = nsin/?/\/2. 

In order to obtain masses of the fermions we diagonalize the mass matrix by unitary L 
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and R matrices: 


LmtR!' = diag(mt^,mi2). 


( 3 ) 


We always set the first eigenvalue equal to the top quark mass, while the second is the mass 
of the new vectorlike quark. 


The mass matrix of the scalars takes the following form: 



( 4 ) 


in the basis $ = where = (^3 — gsin 6 'vi/) cos(2/?)M| is the electroweak D 

term contribution, and A and B are soft breaking terms corresponding to the appropriate 
couplings in the superpotential. Due to mixing with the vectorlike quark, the top Yukawa 
coupling can now be very different from its MSSM value while still keeping the predicted top 
mass unchanged. There are always two values of the top Yukawa that predict the correct 
top mass, and we always chose the larger one. The smaller value is a modihcation of the 
fermiophobic Higgs coupling approach, and generally is more constrained by the data. 

In what follows we consider two sets of new parameters. One set incorporates the small 
mixing example with m = 0 , and the other incorporates the large mixing case with m = 
Msusy- In both cases the superpotential vectorlike mass term M is also equal to Msusy- 
New scalar soft masses are mp = mp = Mgjjgy and all other mass parameters which were 
not present in the MSSM are set to Bm = Bm = Af = 0. For simplicity we set the 
pseudoscalar mass tua and all MSSM soft breaking terms to Msusy except rriHi, and 
B which we vary in order to achieve correct electroweak symmetry breaking for each value 
of Msusy- A-terms are all set to —Msusy- As mentioned above Yt is always hxed by 
requiring that at the tree level rrit^ = which corresponds to the physical top mass 

rrit = 173.35 GeV when one-loop corrections are included. The only free parameters left are 
Msusy and tan f3. 

A. Higgs mass correction 

We calculate the contribution to the mass of the light neutral Higgs boson using effective 
potential approximation in the decoupling regime [5]. The contribution to the effective 
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potential from tops and stops and the new vectorlike states reads 


Al/ = 


647r^ 


i=l 


( 5 ) 


where F{x) = x‘^\n{x/Q'^) while M^, and m? are eigenvalues of the fermion mass matrix ([^ 
and scalar mass matrix (|^ respectively. The correction to the light Higgs boson squared 
mass is equal to 


Ami = 


sin^ d 


dvi 


LA. 

Vn dVn 


+ 


cos^ /3 


Id. . 

^-^ + sin d cos d 

dv^ Vddvd 


52 


dvddv^ 


AH. (6) 


Since the above correction already includes the top and stop contribution, we subtract 
the MSSM top and stop correction which was already included (among other 

corrections lEl) in our MSSM value We calculate the correction using 

eigenvalues of the MSSM mass matrices in equation ([^ and then using an equation similar 
to (§. with only MSSM masses. Our hnal computation of the corrected Higgs mass reads 


mi = (mf + Ami “ (Amf 


MSSM\2 


(7) 


Figure shows the value of Msusy needed to obtain = 125 GeV as a function of 1^/ 
together with various constraints explained in the following section. Figure shows the 
minimal value of Msusy achievable without violating any of the experimental constraints. 
The smaller the value of Msusy the more the vectorlike extension of the MSSM helps to 
ameliorate the little hierarchy problem. The MSSM values of Msusy corresponding to 
tan/3 = 5,7,10 and 30 are Msusy = 11.4,7.4,5.7 and 4.4 TeV, which means that in all 
presented cases we are able to achieve much lower Msusy than required in the MSSM, 
without violating any of the constraints. 

Since the additional contribution to the Higgs mass from the vectorlike quark sector 
lowers the value of Msusy needed to achieve the observed Higgs mass, it also increases the 
prospects of hnding the correspondingly lower superpartner masses at subsequent runs of 
the LHC . 


B. RGE corrections 

The introduction of additional states and additional Yukawa couplings to the MSSM 
causes the renormalization group flow trajectories of the couplings to be altered as the scale 
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FIG. 1: Common superpartner mass MsusY required to obtain rrih = 125 GeV as a function of Yf 
for m = M (left panel) and m = 0 (right panel). Bottom row shows a zoom of the top row plots‘ 
lower right corners. MSSM values of Msusy required to obtain nifi = 125 GeV corresponding 
to tan/3 = 5,7,10 and 30 are Msusy = 11.4,7.4,5.7 and 4.4TeV. Dashed lines are allowed 
by all considered constraints, while solid lines correspond to different exclusions which will be 


achievable in HL-LHC. The calculation of these bounds is explained in section III Dark blue 
regions may be excluded by measurement of the Higgs boson signal strength at 2a significance. 
Dark green regions predict corrections to oblique parameters that may be excluded by future HL- 
LHC measurements at 2a significance, and red regions may be excluded in the second LHC run 
by direct detection of the top partner. Vertical lines show maximal Yu allowing gauge coupling 
unification before the quasifixed point sets in. All parameters except tan /3 are fixed by assuming 
a single supersymmetry scale Msusy and requiring correct top and Higgs physical masses mt = 
173.35 GeV , mu = 125 GeV . 


increases. In this section we discuss these effects and specify the implications and constraints 
they have on the unification of couplings and the possible development of Landau poles in 
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FIG. 2: Minimal value of Msusy achievable without violating any of the above constraints as a 
function of m^sy panel) and tan (5 (right panel). All other parameters are fixed by assuming a 
single supersymmetry scale Msusy and requiring correct top and Higgs masses mt = 173.35 GeV , 
ruh = 125 GeV . 

the couplings. 

In this analysis we have calculated two-loop renormalization group equations using 
SARAH [13], and confirmed the results analytically using known results [T3|. Very sig¬ 
nificant changes in the renormalization group trajectories come from new coefficients in the 
one-loop running of the gauge couplings, 

^ '’■ft = (j. 1. -2) ■ (8) 

These new equations predict the unihcation scale Mu (dehned here by gi{Mu) = g 2 {Mu)) 
to be signihcantly lower than in the MSSM. The new unihcation scale is not far above 
10^3 GeV. 

It is important to point out that unihcation at a scale around 10^® GeV can still easily be 
achieved by positing appropriate high-scale threshold corrections [13] or by adding vectorlike 
quarks so that together all vectorlike superhelds form a complete representation of SU{5). 
This can reestablish coupling constant unihcation without signihcant modihcations to other 
bounds discussed in the following sections. 

However a more stringent constraint comes from the running of and its contribution to 
the running of Y^. At one-loop order these contributions induce Landau poles in the Yukawa 
couplings’ running when V' is sufficiently large — at two-loop order V and Yti develop a 


7 










strongly coupled UV quasifixed point. The range of values of Yti that allow gauge coupling 
unihcation before the UV quasihxed point sets in are Yf G (—1.775,0.002) for m = Msusy 
and Yt! G (—0.8275, 0.8275) for m = 0. These values are marked on the plots showing our 
results. However, since we do not consider a specihc UV completion, it is not necessary to 
treat them as constraints. 

III. CONSTRAINTS 

A. Oblique parameter corrections 

We calculate the S and T parameter |T6] contributions from the vectorlike quarks and 
their scalar superpartners using results from [5], details are shown in Appendix]^ To cal¬ 
culate MSSM contributions we use expressions from Cl excluding corrections from stops 
and sbottoms which were already included in the vectorlike contribution calculation. We 
verihed dominant corrections coming from new fermions with similar results from na. 

The currently allowed experimental values are S = 0.06 ± 0.09 and T = 0.1 ± 0.07 
(assuming U = 0) with correlation 0.91 [19] (the correlation parameter is the tilt in the 
ellipse in the S-T plane). Only minimally more stringent constraints can be achieved from 
LHC running at ^/s = 14TeV with high integrated luminosity 300fb~^. Predicted future 
sensitivity values of S' = 0.06 ± 0.09 and T = 0.1 ± 0.06 are taken from [20] . 

Figure]^ shows resulting corrections to the T parameter as a function of Yf together with 
points showing values above which the results can be excluded at 2 ct by future experimental 
constraints. These points are very close to forming a vertical line because corrections to the 
S parameter are very small for all interesting values of F)/. This is also the reason for which 
we do not include a plot of vectorlike corrections in the S- T plane. 

Corrections from other superpartners are very small due to the simplihed spectrum we 
chose. Figure]^ shows corrections coming from MSSM with and without the stops contribu¬ 
tion from 100,000 randomized spectra of masses up to 3 TeV. A more randomized spectrum 
is unlikely to produce points outside the the S and T exclusion ellipse. Most of the points 
would bring our results closer to the central values due to negative T competing against large 
positive vectorlike quark corrections and positive S contributions, which push our results 
towards the experimentally allowed ellipsis. 




- tan;S=5 

- tan;S=7 

- tanjS=10 
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FIG. 3: Correction to the T parameter as a function of Y^i for m = M (left panel) and m = 0 
(right panel). All values satisfy rrih = 125 GeV . Green points show values above which the results 
can be excluded at 2cj by future experimental constraints. 

Superpartner corrections to electroweak precision observables are generally small because 
superpartners are largely decoupled even with current direct detection exclusions. However 
inclusion of a new quark can introduce unacceptably large corrections to the T parameter 
if its mixing with the SM top is substantial. Nevertheless, it is important to note that with 
currently available bounds, electroweak corrections are the most important constraints on 
our model. However, as the energy and luminosity increase for HL-LHC the observables at 
play in the electroweak precision analysis do not improve substantially. Therefore, precision 
electroweak analysis constraints become relatively less important in time compared to direct 
detection probes of new states and especially compared to precision Higgs analysis, which 
is discussed in the next section. 

B. Higgs boson coupling corrections 

Next we turn to calculation of Higgs boson branching ratios including the above modihca- 
tions and new couplings to the top quark and its vectorlike partner. We start by discussing 
the shifts in couplings of the MSSM compared to the SM and then compare with the case 
with extra vectorlike top states. In the MSSM, the Higgs couplings to up and down type 
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FIG. 4; Oblique parameter corrections in S' —T plane coming from the MSSM (left panel) and the 
same results without stop and sbottom contribution (right panel), with a randomized spectrum of 
superpartner masses up to 3TeV. 


quarks and vector bosons take the form PH 

Qu cos a 

9u^ sin/3 

_9^ ^ -sing . . 

gf^ cos/3 ^ ’ 

9v ■ (o ^ 

9v 

where a is the Higgs mixing angle and tan /3 = Vu/vd- 

Most experimentally important branching ratios have the same values as in the MSSM, 
which are obtained by multiplying the appropriate q coefficients in front of the SM partial 
width exprressoins 

T{h bb) = clT^^ih bb), T{h Tf) = clT^^ih Tf), 

T{h —)■ /i/i) = c^r®'^(h —>■ fip), T{h —)■ cc) = c^r®'^(h —)■ cc), (10) 

T{h ww) = 4r™(h ^ ww), T{h zz) = 4r®^(h ^ zz). 


The remaining important branching ratios are loop induced and are modihed due to modihed 


top couplings and new particles in the loops. We will express these branching ratios as 


V{h^X) 


i^xr 




( 11 ) 


In the following Nc = 3 and loop functions F, I and A, as well as coefficients r, are dehned 
in pT]. Charges and third components of isospin for helds used in the following equations 
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/ 

ti u d e 

e/ 

2 2 1 1 

3 3 3 

Rf 

-‘s 

1 r 1 1 1 

2 nil 2 2 2 


TABLE I: Charges and effective third isospin components. The mixing matrix L is defined in Q. 

are shown in Table |T| while modihcations of the top and top prime conplings to the Higgs 
bosons are given by 


QhtiU 


YtLiiRii + YfLiiR 


i2 


y^MSSM 


( 12 ) 


where L and R are fermion mixing matrices dehned in (|^. are snms of the contri- 

2 

bntions of superpartners which we neglect since they have very small couplings g ~ 

IVl or 

For branching ratio to two gluons we have, 


^SUSY 


SUSY 
99 ’ 


•^99 = Cd Fi{Ti) + Cu Fi(r*) + 9huu^i{Tu) + Al; 

i=d,s,b i=u,c 2=1 

= E E 

i=d,s,b i=u,c,t 

Similarly for the branching ratio to two photons we have, 

Ay-, = cvFi(tw)- t-Cdel ^ FUTi) + cjNcel ^ f i (n) + Fi (n) 


(13) 


i=e,fi,T 


i=d,s,b 


i=u,c 


+ Cu^c^u 'y ^ QhtitiFl ('tJ + 


SUSY 

'77 


(14) 


2=1 


= Fy(Tw)+el Y, F,_{T,) + NY‘d Y Uin)+ Nyel Y UW- 

i=e,iJ,,T i=d,s^b i=u,c,t 

Lastly for branching ratio of Higgs to a photon and Z boson we obtain 


i=e,fi,T 

2 


i=d^s^b 


+ CuNcCu vughtitAh(Fy, AiJ + cvAi{tw, Xw) + 


SUSY 


(15) 


2 = 1 


= CeVe y] Hi(ri,Ai)+ iVce^ Hi (r*. A*) + iV^e^ Hi (r*, Aj) + Hi(riy, A^/), 


i=e,^,T i=d,s,b i=u,c,t 

where Vf = {2T^ — icfS^r)/{swcw), sw = sin 9w and cw = cos 9w- 
The branching ratios are given by 

Tx 


B{h -^X) = 


EC 


(16) 
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<1 


^fJjg-T 

Agww 

Apzz 

0.06 

0.11 

0.08 

0.06 

0.07 


TABLE II: Higgs signal strength future experimental sensitivities at la significance from CMS j23] 


with the sum running over all decay channels computed in this section. We approximate the 
resulting signal strength modification by including only the gluon fusion production channel, 
which at leading order gives 

aB(h ^ X) - X) aB{h X) 


A/ix = 


^SMBSM^h ^ X) 
a{gg —)■ h) Br{h —)■ X) 


- 1 


(TSM5SM(/i ^ X) ^ 

nh^gg) B{h^X) 


(17) 


- 1 . 


a^M(^gg h) B^^{h ^ X) ^ ^ gg) ^ X) 

We confront these results with future experimental bounds as predicted by the CMS 
Collaboration [23] shown in Table SM values of the branching ratios were taken from 
PI. The resulting signal strength modihcations are dominated by the increased gg ^ H 
production cross section compared to the SM and even MSSM. In our model all signal 
strengths grow rapidly when the mixing with the vectorlike state is increased. The most 
important exclusion limit comes from the H —)■ WW signal. The high sensitivity in this 
channel is due to the onset of high statistics and high accuracy in the measurement of this 
channel at the HL-LHC. This can be compared to if —)■ 77 which is not as useful due to 
smaller modifications of its total a ■ B rate. There is a partial cancellation of vectorlike 
top contribution in the a ■ B product. The second best exclusion channel is ii ZZ with 
slightly worse experimental accuracy. The increased experimental sensitivities at HL-LHC 
leads to the conclusion that the hrst evidence for vectorlike quarks in this context of natural 
supersymmetry would likely come from deviations found in precision Higgs observables. 


C. Direct Detection 

The best source for the direct mass bound for the new vectorlike states are dedicated 
analyses by the ATLAS and CMS collaborations at LHC. In particular, the recent CMS 
analysis of t' decaying in three channels t' —)■ bW, tZ, tH without assumptions on the 
branching ratios, has current mass limits between 687 GeV and 782 GeV. 

A similar analysis of decay to the same final states in future colliders was performed 
in |26|. The authors predict mass ranges in which t' could be discovered or excluded for 
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FIG. 5: Vectorlike top partner mass for which mh = 125 GeV as a function of 1^/ for m = M (left 
panel) and m = 0 (right panel). Horizontal line corresponds to the future experimental bound. 

different energies and integrated luminosities. We use their exclusion limit (at 95% C.L.) 
for vectorlike top partner achievable in LHC at ^/s = 14TeV with integrated luminosity 
300 fb“^, namely rrit' < 1525 GeV. 

Figure shows the vectorlike top partner mass needed to achieve rrih = 125 GeV as a 
function of F)/. The right-hand side plot is very similar to Figure because, as expected, 
the mass of the vectorlike top is close to Msusy, while in the left-hand side plot the mass is 
significantly enhanced due to large mixing. 

It is important to point out here that direct detection is crucially dependent on the mass 
of the additional quark, while all previously discussed constraints were more dependent on 
its mixing with already observed states. Gonsequently the interplay between constraints 
described in this section and those of the previous two depends on the mixing, which is 
a consequence of our choice of spectrum parameters. This is why we include both small 
(m = 0) and maximal (m = Msusy) mixing scenarios in our analysis. Direct detection 
bounds turn out to be very important for our model. And in fact this probe proves to be the 
strongest for the part of parameter space corresponding to large tan jS, unless the mixing is 
sufficiently large (m ~ Msusy)- Otherwise precision Higgs analysis will be a more powerful 
probe as shown in Figure 
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FIG. 6: Region of vectorlike mass and mixing parameter space, where direct detection is the 
strongest constraint. The unmarked regions corresponds to precision Higgs measurements being 
the strongest constraint. 

IV. CONCLUSIONS 

In summary we analyzed a single vectorlike top partner model, which is the simplest vec- 
torlike extension of the MSSM that can signihcantly help with the little hierarchy problem. 
We calculated and compared different experimental constraints the model will face after 
300 fb“^ of data are gathered at the HL-LHC. Our key result is that the most constraining 
of the discussed bounds is modihcation of the Higgs boson properties. An exception to that 
is the case of large tan (3 and small mixing where the direct detection probes of the heavy 
vectorlike states at the collider are slightly more stringent. 

After including all the constraints achievable at the HL-LHC, the resulting Msusy can 
still be as low as 1.2 to 2.4 TeV for the simplest possible supersymmetry spectrum. These 
results are 3 to 5 times smaller compared to what otherwise would be allowed in the MSSM. 
Thus even a very simple vectorlike quark extension can greatly reduce the little hierarchy 
problem of the MSSM, and careful measurements of Higgs boson observables would likely 
give hrst evidence of this scenario. 
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Appendix A: Oblique parameter corrections 


The Peskin-Takeuchi precision electroweak parameters [IB] S and T are defined in terms 
of electroweak vector boson self-energies as 

aS 


4s?, -2 


w^w 


C2W 


UzziM^z) - nzz(o) - 


CwSw 

2 TT /n^//\#2 


aT = — Ilzz{0)/Mz 


/Ml (Al) 
(A2) 


The one-loop functions G{x), H{x,y), B{x,y), and F{x,y) have been defined in ref. [T7] . 
Particle names stand for the squared mass of the particles when they appear as an argument 
of these functions. 

Contributions from t' to the electroweak vector boson self-energies are: 




AII^^ — 


Nr 






Nr 


Afl^^ = 


Ibvr^ 


gsw 




i=l,2 


A frSM 


(A3) 


N 


An 


Ibvr^ 

Nr 


ww — 


- 4Re{gf 

i j i j j j 


- An|“, 


tSM 

'-WWi 


i=l 


where Nc = 3, e„ = 2/3, = —1/3 and SM contributions are similar to those above with 

couplings in which Ln = 1 is the only nonzero element of the mixing matrix. The massive 
vector boson couplings with quarks are 


' / Cw \2 / 'J Cw 


w _ _y_T* 

9tib/ 

where L is the fermion mixing matrix defined in ([^. 

The up-type scalar mass matrix (|^ is diagonalized by the unitary matrix U\ 

UMlU^ = diag(m|, m|, m|, m|), 


(A4) 


(AS) 
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while the MSSM sbottom mass matrix is diagonalized by the unitary matrix D 


DM'ljDt = diag(m|^, m|; 


(A6) 


Contributions from third family squarks to the electroweak vector boson self-energies are 

4 2 


Nc 9 ■? 


^UZ'y = 


Nr_ 


167r2 
Nr_ 


gsw 


; Fiii.ii) + el'Y^ F(bi, h 

i=\ i=l 

4 2 

ii) + edY^ 9li*F(hi, hi) 


— 


An 


IGtt^ 

Nr 


ww 


2=1 2=1 

*j=i *j=i 

EEisM-t-fA.t). 


1671-2--.-j 

1=1 j=l 


where the vector boson couplings with the squarks are: 


sli- = “ (y-\(DaDji) - ejsJ,,i 5 yj , 




(A7) 


(AS) 
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